Event-shape selection cuts for supersymmetry searches at the LHC with 7 TeV energy 
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We investigate the prospects of supersymmetry searches at the LHC with 7 TeV energy. A new set 
of selection cuts is proposed based on event shapes to control backgrounds. Our preliminary studies 
show that it is possible to minimize backgrounds to a significantly low level and a conservative 
estimate suggests that the mass reach can be extended to ~1.1 TeV for luminosity £, —IFo with 
a reasonable signal-to-background ratio. 
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INTRODUCTION 



For more than a few decades it has been well known 
that supersymmetry(SUSY) is one of the most promis- 
ing candidates for beyond standard model phsics. Since 
March 2010, LHC has been running with a center-of- 
mass energy 7 TeV and is expected to accumulate lfb~^ 
data by the end of 2011. In both the experiments, AT- 
LAS and CMS are heavily engaged to look for a SUSY 
signal in the collision data. From negative searches 
CMS (ATLAS) have predicted a lower bound on gluino 
mass, nig ~ 650(700)GeV for almost degenerate squark 
mass(mq:) [HQ in the framework of the constrained min- 
imal SUSY model (cMSSM) . It is now of general interest 
to address the discovery potential of SUSY at this low 
energy(7 TeV) and low luminosity(lfb~^) option in the 
LHC experiment. Already quite a few studies have been 
carried out along this direction 0-[5j . 

In this present study, we revisit the discovery reach 
for SUSY at 7 TeV energy with a special effort to control 
standard model(SM) backgrounds by implementing a dif- 
ferent search strategy. In hadron colliders, sparticle pro- 
duction, because of its long cascade decay chain results in 
events with a multiple number of leptons and jets along 
with missing transverse energy(_^y) due to the presence 
of neutralino(xj'), assumed to be the lightest sparticle. 
Since they originate from massive g, q states, therefore 
the final states consist of harder objects. Hence, SUSY is 
characterized by high transverse momentum(pT) events 
with higher multiplicities. We introduce event-shape 
variables @ to exploit these special features of SUSY 
events to discriminate signal and backgrounds. More- 
over, we also try to construct an additional cut based 
on pt of jets in the final state. Feasibility of these new 
sets of cuts are examined by analyzing the following final 
states, which are thought to be discovery channels for 
SUSY at the LHC: 

• a single lepton -I- jets(l^), 

• di-leptons-|- jets(2£), 

• jets + If J,. 

In a hadron collider machine, measurement of Ifrp is a 
nontrivial task due to the presence of other nonphysics 
sources of Ifj,. However, current studies show that Ifj, 
performance in the detector is better than it was thought 
to be 0I- Nevertheless, we investigate the detection pos- 
sibility of SUSY signal with and without I^rp ^ in the 
final state. 





PI 


P2 


P3 


P4 


mo 


500 


1500 


500 


450 




200 


200 


400 


500 


m-g 


524 


575 


954 


1161 


niq 


660 


1535 


981 


1133 


m-± 

^1,2 


142,296 


126,241 


308,515 


391,623 


•^1,2 


78,143 


76,130 


164,309 


207,392 




274,295 


196,240 


499,514 


610,623 


a(pb) 


2.5 


0.32 


0.08 


0.018 



TABLE I: Masses(in GeV) of SUSY particles for four sets of 
mo,mi/2 a-nd fixed values of Ao=0, tan/3=45, sign(/i)=-|-l. 
The total cross sections((T) for SUSY particle production are 
in the last row. 



As is the practice, we simulate signal event in the frame- 
work of a minimal supergravity(mSUGRA) based SUSY 
model described by four parameters, mo, mi/2, tan P 
and sign(/i) at the GUT scale. Here rno:'7ii/2 a-i'^ the 
unified masses of scalars and fermions, respectively, Aq 
is the trilinear coupling, tan /? is the ratio of two vac- 
uum expectation values and /i is the Higgs mass parame- 
ter. The mSUGRA model is very severely constrained by 
many low energy experimental data, like direct bounds on 
sparticles as well as from dark matter experiments [3, ■ 
Instead of testing the mSUGRA model against all those 
constraints, we simply use restricted parameter space ob- 
tained by other authors [4] and select a few sets of pa- 
rameters to a simulate signal. The sparticle masses at 
the electroweak scale are obtained by using the renor- 
malization group evolution performed by SuSpect, and 
decay branching ratios of sparticles are calculated us- 
ing the SUSY-HIT package 0. In Table I, four sets 
of parameters(Pl-P4) are presented for fixed values of 
Ao=0, tan ^=45 and sign(/x)=-|-l. Masses of g and q are 
shown along charginos and neutralinos. In all cases (Pl- 
P4) the lighter chargino state(x]'') turns out to be the 
next-to- lightest supersymmetric particle(NLSP). 



2. SIGNAL AND BACKGROUND 

At the LHC, pairs of colored sparticles, viz gluino- 
gluino, gluino-squark, and squark-squark, are produced 
copiously and the corresponding total leading order (LO) 
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cross sections are presented in the last row in Table I 
for each of the parameter space. We set both the renor- 
malization and factorization scales to = s, defined to 
be the center-of-mass energy in the partonic frame, and 
CTEQ6L [lli] is used for parton distribution function. 
The dominant SM backgrounds are, tt+jets, W/Z + jets, 
ttW+jets, tbW+jets, QCD. In addition, we also check 
contributions due to WW+jets, ZZ+jets and WZ+jets. 
Since, the background to signal cross section ratio are of 
several orders of magnitudes, therefore, to identify the 
SUSY signal, one needs to suppress backgrounds a by 
huge amount, which is not a formidable task by a suit- 
able choice of kinematical selection cuts However, 
a continuous effort is always in process to control SM 
backgrounds to achieve a better signal sensitivity ^T§\. 
For example, a new variable, namely ax, was proposed 
to suppress backgrounds, mainly QCD di-jet events [13]. 
In this current study, our goal is to find a better method 
to deal with SM backgrounds. We propose a new strat- 
egy to suppress backgrounds by using well-known event- 
shape variable, namely, transverse thrust(T). In addition, 
we define another new variable taking the ratio(i?T) of 
a scalar sum of transeverse momentum of lowest number 
of jets over the sum of all jets present in the event. To 
the best of our knowledge these two cuts are never used 
in SUSY searches in the hadron colliders. We observed, 
which are discussed below, that these two cuts play a 
very important role in isolating SM backgrounds leading 
discovery reach signal cross section limited. Importantly, 
these two cuts are also very easy to implement experi- 
mentally once four momenta of jets are reconstructed. 
The event generator PYTHIA6[15| is used to gener- 
ate signal events and background processes due to ti, 
WW,WZ,ZZ and QCD. The ti and QCD backgrounds 
are generated by slicing the entire phase space in various 
Pt bins, where px stands for the transverse momentum of 
final state partons in the partonic center-of-mass frame. 
For QCD, we present results from pr >200 GeV onwards 
since contribution due to low px bins are expected to be 
negligible because of strong cuts as can be seen later. The 
hard scattering process consisting of more than two par- 
ticles in the final state, like tt-|-jets, W/Z-|-iets, ttW-|-jets 
tbW-|-jets, are simulated using ALPGEN fla| based on 
the matrix element (ME) calculation and subsequently 
passed through PYTHIA6 for parton showering(PS). While 
generating events by ALPGEN, initial selection cuts of 
Pt > 20GeV and pseudorapidity, \r/\ < 3 are applied. 
We adopt MLM matching [17] to avoid double counting 
while doing parton showering(PS) after performing ma- 
trix element (ME) calculation. In the MLM matching, 
we use jthe et pr threshold to 20 GeV and |?7| <3 while 
keeping the default value of A_R— 0.7 for jet and parton 
separation. Finally, we multiply matching efficiency with 
the accepted efficiency to obtain event rates. The jets 
are reconstructed by Fast Jet [l^ using the anti-i^T |l9t 
cone algorithm and are preselected with minimum cuts 
p^P >50 GcV and \r]j\ <3. The total I^j, of the event is 
calculated out of the momentum of all visible particles 
present in the event with >1 GeV and {rjl < 3. 



In the following we, describe cuts used in the simulation. 
•Leptons(Cl): Leptons, both electron(e) and muon(/i) 
are selected with pr >W GeV and \ri\ <3. In the case 
of the single lepton final state, we apply pr > 20GeV. 
Isolation of the lepton is ensured by looking at the total 
transverse energy, E^'~' <20% of the pT of the lepton, 
where E^'-^ is the scalar sum of transverse energies of 
jets close to the lepton satisfying AR{i,j) <0.2. 
•Event-Shape variable(C2): These variables describe the 
shape of the event which is related with the geometry of 
the final state and, importantly, are safe to use as ob- 
servables because of its stability against any singularity. 
For instance, the CMS has used event-shape variables to 
test models for QCD multijet production by analyzing 
collision data ^2^]. As discussed before, multiplicity of 
the final state, which determines the shape of the event, 
is very different in SUSY and backgrounds. This encour- 
aging observation leads us to use event-shape variables 
to separate out the signal from the debris of SM back- 
grounds. There are many event-shape variables out of 
which we use only the transverse thrust defined in the 
(x-y) plane transverse to the beam direction, which is 
along the z-axis. The transverse thrust is defined as @, 

T = max„y — =^ , (1) 

l^t iT.i 

where the sum runs over all objects in the event, qT,i 
is the transverse component of each object and nx is 
the transverse vector that maximizes this ratio. No- 
tice that because of the normalization of T by the hard 
scale of the event, the effect of systematic uncertainties 
related with measurements is expected to be less. In 
Fig. 1, in the left panel, we display the distribution for 
T = 1 — T, of signal and the backgrounds for jet plus 
final state. This distribution is subject to preselection 
jet cuts i.e,p^ >50 GeV and \ri\j <3 in conjunction with 
I^rp >50 GeV. It clearly demonstrates that events having 
low multiplicity, like high pT QCD and Z(— >■ i^i/)+iets 
processes, are distributed near t ~ 0(r ^ 1) whereas 
tf-|-jets with higher multiplicities are distributed towards 
comparatively higher(lower) values of r(T). As expected, 
signal events with more multiplicities are distributed to- 
wards larger values of r. This novel feature of transverse 
thrust between signal and background events is exploited 
to suppress the latter. For instance, selecting events by 
demanding T < 0.9(t > 0.1), QCD background and 
W/Z-fjets are suppressed by more than 90% whereas tt 
and the signal are less affected. 

• Rt (C3): Events are accepted requiring a number of 
preselected jets, rij > n™*", where n™*" is the lowest 
number of jets events should have at the least and is the 
input for event selection. Now we construct a ratio(7?T) 
between the scalar sum of pt for the required lowest num- 
ber of jets (n™*") in the event and the scalar sum of pT 
of all jets(nj) in the same event, i.e, 

Rp = Ta_Pt (2) 
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FIG. 1: The transverse thrust(left) and i?T (right) distributions (normahzed to unity) for the signal and the backgrounds with 
selection cuts as shown in the figure. 



with Ht = Y^i' Pt ■ It is obvious that the behavior of 
this variable is strongly connected with the multiplicity 
as well as the hardness of the jets in the event. The Rt 
becomes exact unity for those events with nj = n™™, 
otherwise it is expected to be away from unity for those 
events with rij >> n™*". In Fig. 1, we present in the 
right panel the distribution of Rt of signal along with 
backgrounds subject to jet selection cuts along with 
the requirement of n™™ =4. The peak suggests that 
most of the events in the backgrounds are dominantly 
4-jet events. The dip indicates events with relatively 
softer jets, pip ~50 GeV(recall that preselection cut on 
jets Pt >50 GeV) along with Uj just above n™*", are 
present predominantly in backgrounds). On the other 
hand, for the case nj >> n™™ and for comparatively 
harder jets, which is the case for signal events, Rt 
is mostly distributed below 0.85. Evidently, the Rt 
distribution provides us another very robust tool towards 
the suppression of SM backgrounds by a substantial 
amount. For instance, restricting Rt <0.85 enables to 
remove backgrounds by ~70-80% or or more with a 
10-20% loss in signal events. 

• Ht{GA): As mentioned before, the signal events are 
expected to be harder in comparison to background 
events. Therefore, a cut on Ht >900 GeV is very 
effective to get rid of a good fraction of background 
events. 

• I^rp cut(C5): In all backgrounds, the lepton mainly 
comes from W decay except for QCD, hence the trans- 
verse mass between the lepton and I^j, is expected to 
be bounded by the W mass, for the single lepton case, 
in particular. Hence, events with a single lepton case 
is expected to suffer due to a cut on transverse mass, 

Mt = ^2E^^p{l - cos(t){e, ^p)) >60 GeV (2l|, where 
(j) is the azimuthal angle between the lepton and the 



I^rp direction. Beside this Mt cut, which is applicable 
only for one lepton final state, a cut on I^j, >150 GeV 
for all cases turns out to be a good criterion to reject 
backgrounds. 

• One lepton -|- jets(l^): The event with the final 
state having only a single lepton(e or /i) along with at 
least 3 jets(n™"'=3) are simulated applying a sequential 
set of cuts C1-C5 and the results are presented in Table 
II using PI parameter space for signal. The second 
and third columns display the raw production cross 
sections and the simulated number of events(AfEv)j 
respectively, for each process. The numbers in the 
third column corresponding to tt-|-jets and W/Z-|-jets 
represent the number of events with jet-parton MLM 
matching. We observe that, in general, about 20%-25% 
of the events are lost due to the isolation requirement of 
the lepton, with the exception of the QCD background 
with more(^50%) loss. In the 5th column, the numbers 
after thrust(Eq[l|) cut(C2) T <0.90, clearly indicate that 
the QCD background is eliminated by a huge amount 
whereas other backgrounds suffer by about 50%-70% 
with a modest loss in the signal event. Requirement 
of at least 3 jets with a restriction on Rt < 0.85(C3) 
results in a substantial suppression of backgrounds, 
particularly processes with low jet multiplicity. Sub- 
sequently, a cut(C4) on Ht >900 GeV removes a 
significant fraction of backgrounds. The resulting cross 
sections (cTg^ J,) due to cuts C1-C4 i.e without the I^j, 
requirement turn out to be 179fb for the signal and 
lOlfb for the total background leading to S/VB ~ 17. 
However, for other sets of SUSY parameters P3-P4, this 
ratio is ~ 1. We check that the contributions due to 
backgrounds ttW-|-jets, tbw+jets and VV+jcts(V=W,Z) 
are negligible after applying all cuts. Finally, I^j, is 
considered to suppress the remaining backgrounds to 
gain in signal-to-background ratio. Demanding(C5) 
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TABLE II: Event summary for the signal and the back- 
grounds after each set of cuts(Cl-C5) described in the text 
for the single lepton(l^) case. For QCD and ti+]ets, events 
are simulated for different pr bins as shown. 



the transverse mass(MT) between the lepton and I^rp, 
Mt >60 GeV and >150 GeV it is possible to 
further isolate backgrounds to an almost negligible level ( 
cr^^ ^2.65 fb). Note that due to the cuts C4 and C5, 
the signal suffers substantially because of comparatively 
lighter g and q masses. However, for higher masses, 
the losses are minimal. In Table III, we show the total 
signal(Pl-P4) and background cross sections after cuts 
Cl-C4((To9^) and C1-C5((j^^). This table predicts that 
the single lepton channel alone can explore g and q of 
masses up to ~1 TeV for C =lfb~^ requiring a signal(S) 
to a background(B) ratio, S/Vb > 5. 
• Di-lepton +jets(2^): The final state with two 
Icptons of any type and of any charge combinations 
are accepted along with at least 3 jets(=n™™). In the 
simulation, a softer cut on the lepton, >10 GeV 
(£=e,/j,) is imposed along with the isolation requirement. 
In this case one of the most dominant backgrounds 
appears to be due to Z-|-jets. A good fraction of Z-|-jets 
events can be removed by an additional requirement 
on the di-lepton invariant mass, m^+i- >10 GeV and 
rrii+i^ ^70-120 GeV. We find that the patterns of 
background suppression due to cuts C2-C4 are more or 
less the same as a single lepton case. In Table III, we 
present the total cross section for both the signal and 



the backgrounds with C1-C4 cuts (ctq^^) and including 
l^j, cut(Cl-C5) (cr^^). This table depicts that, without 
Ifrp cut it is difficult to observe any signal event in the 
di-lepton channel except for light g and q masses(Pl). 
Therefore, including I^j, cut(C5) it may be possible to 
detect signal events for g and q masses ~700 GeV with 
£= Ifb-i. 

• Jets +I^rp: It is well known that SUSY searches in 
jets plus the missing energy channel offer the highest 
reach of gluino and squark masses P, We revisit 
this final state to examine the effects of transverse 
thrust (EqllJ and i?T (eql2]) selection cuts in tandem 
with the other set of cuts. A substantial number of 
background events(~50%-80% or more) are rejected by 
thrust cut(C2) costing about 20% of the signal events. 
Furthermore, events with at least 4(=n™") jets with a 
softer cut on the 4th subleading jet, >70 GeV are 
selected and a requirement on Rt <0.85 turns out to 
be very useful to get rid of a good amount of W/Z-|-jets 
and ttW+jets events. However, C1-C5 cuts are not 
enough to kill backgrounds completely, a good amount 
of <t-|-jets remains. Therefore, in addition, we adopt a 
selection cut on the transverse mass to^^"'^ , between two 
leading jets. The demand of mi^^^ >450 GeV is very 
effective in suppressing tt-t-jets events by almost a factor 
10, of course, with a substantial loss('-^60%-70%) of 
signal cross section as well. However, for higher masses 
of g and q i.e for parameter sets P3 and P4, we find this 
loss is less than 5%. The cut Ht >1050 GeV removes 
the remaining background events from W/Z-fjets, 
whereas high px QCD and tt events are eliminated by 
the I^j, >150 GeV cut. We find after all cuts the total 
background cross section mainly due to the ti and QCD 
is ~ 3.7fb out of which tt-fjets contributes about 70%. 
In Table III the total signal (P1-P4) and background 
cross sections are presented. It reflects that at g and 
q masses up to TeV can be explored with lfb~^ 

luminosity. We find that the signal efficiency for low 
mass point (PI) is about 10% and for higher mass(P4) 
it turns out to be --30%-40%. 



3. CONCLUSIONS 

We investigate the detection possibility of the SUSY 
signal in three different channels by implementing a set 
of cuts few of which were not used in earlier analysis. 
Our study shows that one of the eventshape variables, 
namely, transverse thrust, and another observable Rt 
play a very significant role in disentangling backgrounds 
from the signal events. Interestingly, both these ob- 
servables, T and Rt are useful in isolating backgrounds 
consisting of higher multiplicity, unlike the case the 
for aT 14] variable which works for di-jet events only. 
It implies that these variables, T and Rt are compli- 
mentary to aT- The added advantage is that being 
dimensionless variables, measurements of T and Rt 
involve less systematic uncertainties. Our conservative 
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TABLE III: Total signal(Pl-P4) and background cross sec- 
tions(fb) after cuts Cl-C4(cro_9y ) and C1-C5((t^^) for the sin- 
gle lepton(l^), di-lepton(2f) and jets plus l^j, case. 



estimate predicts that in the single lepton channel and as 
well as the jet plus I^j. channel, it is possible to achieve 
a reasonable signal-to-background ratio for g and q 
masses up to ^ 1.1 TeV whereas the di- lepton channel 
alone is not very encouraging. It is to be noted that our 
conclusion is based on LO signal cross sections whereas 
in background evaluation the higher order effects are 
taken into account to a certain extent by considering 
hard emission of partons(jets), which is the real part 
of the NLO correction. We check using Prospino 
that NLO cross sections for SUSY are about 50%-70% 
higher than LO cross sections resulting in an enhance- 



ment in signal-to-background ratio than our present 
conservative estimate. Clearly, the discovery reach is 
signal rate limited rather than background limited, 
which is almost negligible after all cuts. Evidently, this 
is a very powerful result based on our search strategy, 
which can be implemented experimentally very easily. 
For illustration purposes, we mention that the total SM 
background cross section corresponding to the jet plus 
missing energy final state is 3.7 fb, as shown in Table III, 
and which is 2.42 pb for same final states as reported 
in the paper of Ref. [2^. We observe that our analysis 
predicts better sensitivity than others in the litera- 
ture [1, 01 ■ Note that in our simulation t leptons are not 
considered in final states, which will be simulated in the 
future along with a detailed scan of parameter space [131 ■ 
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